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Abstract 

In this paper, a direct adaptive radial basis function (RBF) 
neural network control algorithm is presented for a class of 
ship course with uncertain discrete-time nonlinear systems. 
To aviod some system states that are unmeasurable and 
make the adaptive control approach more universal and 
convenient to be implemented in practical application, the 
original ship course with uncertain discrete-time nonlinear 
system is transformed into the form of the input-output 
model. According to the input-output model, a direct 
adaptive RBF NN control for the ship course with discrete- 
time nonlinear system is carried out based on the existence 
of the implicit desired feedback control (IDFC). In the 
controller design process, RBF neural networks are used to 
emulate the desired desired feedback control and 
approximate the unknown function. The stability of the 
closed-loop system is proven to be uniformly ultimately 
bounded (UUB) by using lyapunov theorem, and tracking 
error can converge to a small neighborhood of zero by 
choosing the design parameters appropriately. In the end, 
the simulation example of motor vessel "yukun" is 
employed to illustrate the effectiveness of the proposed 
algorithm. 
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Introduction 

In the past decades, control over ship motion has 
attracted an ever increasing interest; especially for the 
ship course control, since its properties have a great 
effect on the security and economy of ship navigation. 
Due to the characteristics of ship motion such as large 
inertia, time delay and nonlinear etc., the ship course 
control system has always been a uncertain nonlinear 
system being a complex control problem. In recent 
years, many control methods have been utilized in 
ship course control systems, for instance robust 



1998), backstepping design (Wang, 2002) and adaptive 
control (Liu, 2012). These control schemes are in 
continuous-time domain. 

However, along with the rapid expansion of martime 
industry, the feature requirements of the ship 
maneuverability and control performance is more and 
more demanding. At the same time, ship is becoming 
large-scale, high speed and automation, therefore the 
traditional control scheme could't satisfy the reality 
requires any more. Thus, it is very important that the 
advanced controller for ship course should be studied. 
It is common knowledge that discrete-time systems 
rather than the continuous-time systems are the closest 
for describing a real plant. Following the 
developedment and modification of information 
processing technology and computer technology in the 
modern ship, the control of ship course with nonlinear 
system in discrete-time form has been a expected field 
(Song, 2007; Bai, 2009). 

Since the universal approximation of the neural 
networks (NNs) has been proven (Park, 1991), the 
NNs have become an active research topic and 
obtained widespread attention in the control domain. 
Owing to the excellent function approximation ability , 
NNs have emerged as an attractive and powerful tool 
for stability analysis of complex nonlinear dynamic 
systems (Li, 2004; Ge, 1998; Lewis, 1999). NN control 
has been extensively studied for both continuous-time 
and discrete-time systems (Ge, 2002; Hovakimyan 
2002; Jagannathan, 2006; Ge, 2003). 

In comparison with continuous-time control design, 
NN control design of discrete-time systems is more 
difficult due to the lack of mathematical tools in 
discrete time. For instance, lyapunov design for 
nonlinear discrete-time systems has become much 
more intractable than that for continuous-time systems 
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because the linearity property of the derivative of a 
lyapunov function in continuous time is not present in 
the difference of lyapunov function in discrete time 
(Yeh, 1995; Song, 1993). In recent year, there have been 
many control schemes to be proposed for discrete-time 
nonlinear systems (Yang, 2010; Ge, 2008; Zhang, 2005). 
For example, for a class of a single-input-single-output 
(SISO) affine NARMA system, direct adaptive NN 
control was proposed in (Chen, 1995), then based on 
the implicit function theorem, direct adaptive RBF NN 
control scheme was studied for a class of discrete-time 
single-input single-output non-affine nonlinear 
systems (Zhang, 2002). Furthermore, it was developed 
with multilayer neural network (MNN) for general 
nonlinear autoregressive moving average with 
exogenous inputs (NARMAX) systems in (Ge, 2004). 

Motivated by the previous works, in this paper, a 
direct adaptive RBF NN control would be studied for 
a class of ship course with uncertain discrete-time 
nonlinear system. In order to avoid the problem of 
some system states that are unmeasurable, the original 
ship course with uncertain discrete-time nonlinear 
system is transformed to input-output model. 
According to the proved existing implicit desired 
feedback control (IDFC), the direct adaptive RBF NN 
control scheme has been proposed for ship course 
discrete-time uncertain nonlinear system. The RBF 
neural network are used as emulator of the Implicit 
Desired Feedback Control and approximate unknown 
function. It is proven that all the signals of the closed- 
loop system are uniformly ultimately bounded and the 
tracking errors converge to a bounded compact set. 
Finally, simulation example of the ship "yukun" 
illustrates the effectiveness of the proposed approach 
with desirable preformance. 

Problem Formulation 

System Description 

In general, the relationship between course <p and 

rudder angle 8 in the mathematical model of ship 
course control system can be described as follows (Jia, 
1997): 



1 K 

<p + — H(<p) = — 8 
T T 



(1) 



where K is a constant for rotary of ship, T is a 
constant for trackability of ship. H((p) is a nonlinear 

function on <p , which can be indicated as follow: 



where a i ,i = 1,2,... are positive constants, which 
denotes nonlinear coefficients of the ship. 

Based on formula (1) and (2), the following ship course 
with discrete-time nonlinear systems is obtained 
(Zhang, 2009): 

x^(k + 1) = x 2 (k) 

1 3 K 
x 2 (k + 1) = (a l x 2 (x) + a 2 x 2 (x)) H u(k) (3) 

y(k) = x x {k) 

where x { ,x 2 e R , y(k) e R and u(k) e R are the state 
variables, system output and input respectively. 

According to literature (Liu, 2013), this discrete-time 
system (3) can be described as the following form: 

y(k + \) = f(y(k)Mk)) (4) 
this model with relation of an input sequence 
{«(&)} to an output sequence |y(&)} by nonlinear 

difference equation, constitutes an extremely broad 
class, including many other classes of nonlinear 
discrete-time model as special cases. 

The control objective is to design a direct RBF NN 
controller for system (3) such that the system output 

follows the desired reference signal y d (k) . In the 

following , the following assumptions have been made 
based on the systems (3). 

Assumption 1 The nonlinear function f (/) is 
unknown, continuous and differentiable. 

Assumption 2 The neural number of hidden neurons 
is /. 

Assumption 3 Assume that the partial derivative 

du 

constants 



> s > , where both £ and g± are positive 



This assumption states that the partial derivative 



du 



3 5 
H \<p) = a x (p + a 2 (p +a 3 4> +• 



(2) 



is either positive or negative. From now onwards, 

df 

without lose of generality, it is supposed that — > . 

du 

Lemma 1 (Mean Value Theorem). If a function f (x) 
is continuous on the closed internal [a, b] and 
differentiable on the open interval (a, b), then a point c 
exists in (a, b) such that: 
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f(b) = f(a) + (b-c)f'(c)\ c ^ 



b) 



(5) 



Radial Basis Function NN 



In control engineering, RBF neural networks ususlly 
are used as a tool to model nonlinear functions 
because of their good capabilities in function 
approximation. The RBF can be considered as a two- 
layer network in which the hidden layer performs a 
fixed nonlinear transformation with no adjustable 
parameters, i.e., the input space is mapped into a new 
space. The output layer then combines the outputs in 
the latter space linearly. Therefore, they belong to a 
class of linearly parameterized networks. In this paper, 
the following RBF is used to approximate the 

continuous function s(z) '■ R q —> R : 



S nr,(z) = W h(z) 



(6) 



where the input vector z s D. z a R q , where q is the 
neural network input dimension. Weight vector 
w = [w± , w 2 , • • • , w, ] e R , the NN node number / > 1 ; 

T 

and h(z) = [h l (z),---,h l (z)] , with h-{z) chosen as the 
commonly used Gaussian function, which is in the 
following form 



hj(z) = exp 



-(z-Mi) (Z~Mi) 
2 



,i = 1,2,...,/ 



(7) 



where ju t =[ju n ,ju i2 ,---,ju iq ] is the center of the 

receptive field and J] i is the width of the Gaussian 
function. 

It has been proven that network (6) can approximate 
any continuous function over a compact set Q z 
to arbitrary accuracy as 



*T 



s(z) = w h{z) + s 7 ,\/z efl 7 



(8) 



where w is ideal constant weights, and £ z is the 
approximation error. 

* 

The ideal weight vector w is an "artificial" quan-tity 
required for analytical purposes, w is defined as the 

value of w that minimizes \s z \ for all Z £ ^ z in a 
compact region, i.e. 



w := arg min < sup 



weR 



s{z)-w h(z) 



(9) 



Controller Design and Stability Analysis 

Assume that y d (k + 1) is the system's desired output 



at time k + 1 . In the ideal case, there is no disturbance; 
and it can be shown that if the control input u*(k) 
satisfies 

f(y(k),u*(k))-y d (k + l) = (10) 

then the system's output tracking error will converge 
to zero. 

* 

Definition 1 If a controller u (k) satisfying equation 
(10) exists, then the controller will drive the system 
output to the desired output, in addition, control input 

u*(k) is called Implicit Desired Feedback Control 
(IDFC). 

the tracking error is defined as 

e(k) = y(k)-y d (k) (11) 

and then the tracking error dynamic equation is given 
as follows: 

e(k+l) = f(y(k)Mk))-y d (k+l) 



(12) 



For the system (4), a RBF neural network was 
designed to realize the controller directly [24]. The 

ideal control input is u , and 

u\k) = u\z) = w* T h(z) + s u (z) (13) 
on the compact set D. z , the ideal neural network 

weights w and the approximation error are bounded 
by 



w 



< w 



(14) 



where w m and S\ positive constants. 

w(k) is defined as the actual neural network weight 
value, and the control law is designed by using RBF 
neural network directly 

u(k) = w T (k)h(z) (15) 
where h(z) is Gaussian function output and z(k) is 
input of RBF. 

By noticing (13), we have 

u(k) -u(k) = w T (k)h(z) - (w* T h(z) + e u (z)) 

= w T (k)h(z)-s u (z) 

~ * 

where w(k) = w(k) —w is the weight approximation 
error. 

The weight update law was designed as follows [24]: 

w(k + 1) = w(k) - y{h{z)e{k + 1) + aw(k)) (17) 
where ^ > is a positive number and <T > . 



(16) 
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Subtracting w to both sides of (16), we have 

w(k + 1) = w(k) - y(h(z)e(k + 1) + aw(k)) (is) 

Using mean value theorem, let a = u (k) , C = g , 
b = u(k) = u (k)+w (k)h(z) - £ u (z) . noticing Eqs. 
(10) and (16), then 

f(y(k),u(k)) = f(y(k),u(k) + w T (k)h(z)-e u (z)) 



y d {k + \) + {w (k)h(z)-s u (z))f u 



where f u 



du 



,g <e(u (k),u(k)) . 



(19) 



(20) 



then we get 

e(k + l) = f(y(k)Mk))-y d (k + V) 

= (w T (k)h( Z )-£ u (z))f u 

T e(k + \) 
and w (k)h(z) = + s {z) 

fu 

The stability analysis of the closed system is given as 
follows: 

Choose the Lyapunov function candidate 

V(k) = —e 2 (k) + -w(k) T w(k) 
§l 7 

then we have 
AV=V(k+l)-V(k) 

= — (e 2 (k + 1) - e\k)) --w(k) T w(k) 

Si r 

+ -w(k+l) T w(k + l) 

r 

= —(e 2 (k + l)-e 2 (k))--w(k) T w(k) 



(21) 



S, 



7 



1 



+ -(w(k) - 7(h(z)e(k + 1) + ow(k))) 

7 

x (w(k) - 7(h(z)e(k + 1) + ow{k))) 
= —(e 2 (k + 1) - e 2 (k) - 2w(k) T h(z)e(k + 1) 



S, 



2aw(ky w(k) + 7I1 1 (z)h(z)e z (k + 1) 



+ lyow (k)h(z)e(k + 1) + ya L w{k) L w(k) 



2 \T ~ . 



Since 



h T (z)h(z)= h(z) 2 </,/ = !, 2,...,/ 



vv 



) 



2avv(£) T = cr(||iv(fc)|| + ||w(*)|f 

^ r (z)h(z)e 2 (k + 1) < ^/e 2 (Jfc + 1) 
2j/aw r (k)h(z)e(k + 1) < yerf[|w(*)|| 2 + e 2 (Jfc + 1)] 
;Kcr 2 w r (fc)vv(fc) = ^cr 2 ||w(fc)|| 2 



We obtain 
Si 



AV < — (e 2 (k + 1) - e 2 (fc)) + ?/o- 2 to)|| 2 



_ 2( ^+l) + ^ (z))e(fc + 1) 

2 2 2 

- cr(||w(fc)|| + - w* ) 

+ jie 2 (k + 1) + ^[||w(£)| 2 + e 2 (fc + 1)] 

< (— - A + r (l + a)l)e 2 (k + l)- — e 2 (k) 
Si Ju Si 

II II 2 II *ll 2 

-2s u (z)e(k + 1) - a\\w(k)\\ +<j\\w 

+ cr(-3 + 7I + 7<j)\\w(k)f 

Based on the assumption (13), from < e < f u < g l , it 
can be deduced that 

Si f u Si Si Si 
- 2s „ (z)e(k + l)< k Q sf + — e 2 (k + 1) 



where k Q is a positive number. Thus, we have 

Si K 



AV < ( — — + /(l + <j)l + —)e 2 (k + 1) 



2 1 

+ cr(^ + 7cr-l)||w'(A:)|| e 2 (fc) 

Si 

- cr||w(A:)|| 2 + ow 2 n + k Q sf 
= -{—-{\ + cj)ly-±-)e 2 {k+\) 

Si K 

+ a((l + <r)y - 3)|w(£)f - a\\w(k)f 

-(e 2 (k)-j3)/ gi 
where (3 is a positive number as 

P = Si(^+Vi 2 ) 



(22) 



Choosing the positive constants as k Q , >l and CT , 
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these contants must be satisfied the following 
inequalities: 

— - — >0 Q. + o)ly >0 

.e, k n g, k n 



(l + a)y-3<0 



that is, 



< g 1 < k Q 

1 1 

Si k o 
0<(l + a)y<3 



0<(l + a)ly< 



(23) 

(24) 
(25) 



then we have AV < once c (&)>/?. This states 

that for all k > , V(k) is bounded because 

k 

V(k) = V(0) + Z AV(i) < oo 
7=0 

Define compact set Q e = |e| e 2 < yffj , then it can be 
seen that the tracking error e(k) will converge to Q e 
if ^(k) is out of compact Q e . 

Simulation Study 

To illustrate the effectiveness of the proposed control 
scheme, the motor vessel "yukun" of Dalian Maritime 
University is used to implement the simulation. 

Consider the ship course discrete-time nonlinear 
system as: 

x { (k +1) = x 2 (k) 

1 3 K 

x 2 (k + 1) = - — {a l x 1 {x) + a 2 x 2 {x)) + — u(k) 
y(k) = x 1 (k) 

the ship' length L is 116 m, breadth B is 18 m, draft 

d is 5.4 m, and the speed is 16.7 kn. Therefore, the 
mathematical model parameters of ship's nonlinear 
motion is K =0.5 , T = 64 , a x = 0.5 , a 2 = 30 . The 
tracking objective is to make the output y follow a 
desired reference signal that is y m (k) = 30sin(£;r/l80) . 

Since g 1 > df fdu = l/64, thus we can set g x = 0.02, 
and from (23), we can choose k^ = 0.05 . For RBF 
neural network, the NN node number / = 2 1 . The 
initial value of the plant is x(0) = [0, 0] r . To satisfy (24) 
and (25), we can choose y = 0. 1 1 , a = 0.9 . 



The simulation results are shown in Figs. 1, 2 and 3. 
Fig. 1 shows the tracking performance for output y 

relative to the reference signal y^ik) obtained. Fig. 2 

and Fig. 3 illustrate the control input signal of the 
closed-loop system and the norm of RBF neural 
network weight respectively. 
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Conclusion 

In this paper, a class of direct adaptive RBF neural 
network controller is presented for the ship course 
with uncertain discrete-time nonlinear system. 
Motivated by the previous works, the ship course with 
discrete-time nonlinear system is transformed to 
input-output model, in which the existance and 
uniqueness of implicit desired feedback control have 
been proven by Implicit function theorem. Then based 
on the input-output model, RBF neural networks are 
used as emulator of the Implicit Desired Feedback 
Control. All the signals of the closed-loop system are 
guaranteed to be uniformly uniformly ultimately 
bounded, and the tracking error can be made 
arbitrarily small by choosing the parameters in the 
control law. Simulation example of the real ship 
"yukun" has been carried out to demonstrate the 
performance and the effectiveness of the proposed 
algorithm. 
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